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Abstract: Following the nonperturbative prescription for the jet quenching parameter recently 
proposed by Liu, Rajagopal and Wiedemann, we compute the first correction in the inverse 't 
Hooft coupling corresponding to string a' corrections in the dual background. We also consider the 
introduction of a chemical potential for the U(l) 3 gauged R-symmetry. While the former mildly 
diminishes the jet quenching parameter -this suggesting a smooth interpolation between the strong 
coupling and perturbative results-, the latter generically increases its value. We comment on the 
extension of this setup to quarks of finite mass. 
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1. Introduction and Summary 

The experimental program at the Relativistic Heavy Ion Collider (RHIC) in Brookhaven National 
Laboratory |I[], has provided considerable insight into the properties of strongly interacting matter 
at high energy densities. Phenomenological analysis has established several striking features of such 
substance. First, the results on elliptic flow are well described by hydrodynamical models only if 
the shear viscosity is taken very low. The medium behaves like a strongly coupled plasma which 
resembles a liquid more than the gas of quasi-free partons long expected to be the state of matter 
at such energy densities, known as the Quark Gluon Plasma (QGP). As a second indication of this 
behavior, high energy partons traversing the medium are strongly quenched. This phenomenon is 
usually characterized, in models of medium-induced radiation, by the so-called quenching parameter 
(or transport coefficient) q @. This parameter has the meaning of the average squared transverse 
momentum transferred from the medium to the traversing parton, per unit mean free path (see the 
reviews 0, [J). 

Phenomenological models differ in the detailed framework for calculating the radiative energy 
loss 0, ||, in the treatment of the geometry and dilution of the medium ||, as well as in the 
consideration of flow- induced radiation and of additional elastic scattering ||. The extracted 
values of the transport coefficient are q ~ 1 -f- 15 GeV 2 /fm, substantially larger than those found 
in studies of hadron production in DIS on nuclear targets, see e.g. |§. While the lower bound 
is compatible with expectations from perturbative QCD [[It]], higher values demand additional 
non-perturbative mechanisms. Therefore, and while waiting for upcoming both experimental and 
phenomenological efforts, it is of uttermost importance to get further information on the possible 
values of q in the strong coupling limit. 

The traditional tool for such studies, namely lattice QCD, cannot be presently applied to 
determining the jet quenching parameter. In contrast, AdS/CFT duality |TT| provides a power- 
ful calculational framework where quantum properties of supersymmetric Yang-Mills theories at 
strong 't Hooft coupling A = g\ M N c and large number of colors N c , are translated into classical 
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computations in a gravitational background. The applications of AdS/CFT techniques to thermal 
deformations of gauge theories started in [I5|. There, the thermodynamics of the black brane ge- 
ometry was conjectured to describe the behavior of the dual quantum field theory at the Hawking 
temperature of the black hole. In the limit of a flat horizon, the dual field theory lives in an un- 
confined phase at strong coupling. Placing probe charges in such a background corresponds to the 
insertion of sources for very massive quarks in the dual theory. Following this lore, a number of 
important results where derived concerning the qq potential, including features like confinement and 
screening both at zero and finite temperature ]n|. They typically involve a Wilson line stretching 
either in a timelike or a spacelike direction; rotating Wilson lines where also examined as putative 
duals to high spin mesons. More recently, the study of thermodynamical properties was extended to 
encompass near equilibrium magnitudes [14]. A not minor surprise came out with the finding of a 
universal ratio between the shear viscosity and the entropy density, r]/s = 1/ (47r) |L5| for quantum 
field theories admitting a holographic dual description. This ratio was conjectured to set a universal 
bound on physical thermal field theories. The data at RHIC suggest that the values for the QGP 
are compatible with the lower bound, this strongly supporting the use of AdS/CFT to describe 
such a system. 

Motivated by these successful applications of AdS/CFT to the study of strongly coupled phe- 
nomena in thermal gauge theories, Liu, Rajagopal and Wiedemann (LRW) recently proposed a 
scheme to determine the jet quenching parameter HTjJ 1 . In their construction, it is central to use 
the identification M of this parameter with the coefficient in the exponent of an adjoint Wilson 
loop computed along a rectangular contour C with a large distance L~ along the light-cone, and a 
spacelike separation L, L <^ L~: 



{W A {C)) = exp 



At large N c this Wilson loop can be expressed in terms of the Wilson loop for the fundamental 
represention, (W A (C)) ~ (W F (C)) 2 . In turn, the AdS/CFT correspondence tells us that this 
fundamental Wilson loop can be computed [O] evaluating the classical Nambu-Goto action S for 



a string ending on the boundary along the previous contour, 

(W F (C)) = exp [S(C)} 



;i.2) 



The result in [IE|, obtained in a near extremal D3 background corresponding to M = 4 SUSY QCD 
at finite temperature, exhibited some interesting features. 2 The quenching parameter q turned out 
to be proportional to T 3 (which of course provides the correct dimensions) and to vA (thus to \/N^), 

1 After the initial proposal, a host of papers have appeared jl7| |l8|, [H]. ||(| |2^] which have shown that this result 
is not universal. Also the very interesting and related problem of a drag force on the brane has been addressed by 
several authors [^2j . 

2 Although the relation between the results in this framework and real QCD is unclear, the hope is that an 
understanding of the size and dependences of the jet quenching parameter could provide some kind of upper bound, 
while the perturbative QCD results should provide a lower bound. In this way, computing the results in less 
supersymmetric backgrounds might provide an indication (assuming a smooth behavior for the transition from 
Af = 4 SUSY QCD to real QCD) of the expected value for this parameter at strong coupling. 
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the latter being totally different to the a priori expected dependence on the number of degrees of 
freedom of the energy or entropy densities, hence oc N%. In this way, the quenching parameter 
appears not to be a direct measure of the energy density of the system, q oc e 3//4 as usually assumed 
TCfl , but of the third power of the temperature. Moreover, the numerical values turned out to be 
astonishingly close to the experimental data: for standard values a s = ~ and N c = 3, q = 0.94 (3.16) 
GeV 2 /fm for T = 200 (300) MeV. Notice, however, that these quantities imply a 't Hooft coupling 
A = 6tt, while the gravity computation is strictly valid in the limit A — > oo. Thus, a more precise 
theoretical value of q in the strong coupling limit neatly demands the understanding of finite 't 
Hooft coupling corrections to the result in JIB] . This is one of the targets in this letter. 

The paper is organized as follows. In the next section, trying to keep ourselves as generic as 
possible, we provide a formula for the jet quenching parameter that can be readily applied to a large 
class of metrics. As a quick example, we present the results for the thermal deformation of Witten's 
D4-brane background Jl2| . We also comment on the fact that this formula admits a straightforward 
generalization to encompass quark sources of finite mass, and provide some preliminar numerical 
analysis. In section 3 we compute the first correction in the inverse 't Hooft coupling to the value 
given by LRW. We show that this correction mildly diminishes the jet quenching parameter. This 
suggests a smooth interpolation between the strong coupling regime and the perturbative results, in 
analogy with what has been observed for the free energy and the ratio rj/s P3], |2~3}|. In section 4, the 
effect of turning on chemical potentials 3 is thoroughly investigated. The relevant metric corresponds 
to the background of a stack of rotating D3-branes with maximal number of angular momenta. We 
explore the evolution of the jet quenching parameter along the space of these three independent 
charges. We typically find an enhancement within the range of thermodynamical stability. We 



further compare with recent results which have appeared on the subject [19, EO, pi 



2. The jet quenching parameter 

In this section we shall provide a formula that allows to readily compute q in string theory back- 
grounds within the class of metrics that are suitable for an AdS/CFT duality, including the case in 
LRW. We will follow essentially the same steps as in [|T7| . The family of ten dimensional metrics of 
interest for us adopt the following form: 

ds 2 = G mn dX^ dX^ 

= -c 2 T dt 2 + c 2 x dx l dxi + 4 dr 2 + G Mn dX M dX n , (2.1) 

where X M = (t, x l , r; X n ), i — 1, ...,p, n = 1, 8 — p. This class of metrics encompasses rotating 
backgrounds which we shall analyze later. We are interested in black brane solutions. 
Following [[16] we will consider the following lightlike Wilson line 

x~ = t , x 2 = a , r = r(a) , (2.2) 

3 These chemical potentials are conjugated to i?-charge densities of Af — 4 SYM theory. They should not be 
confused with that corresponding to the baryon density in QCD whose implementation in the dual supergravity 
side is currently an open problem. Indeed, the baryonic charge is not dual to a U(l) isometry of the supergravity 
background. We thank Krishna Rajagopal for stressing the importance of this point. 
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with x ± = (x° ± x r )/^/2, r G (0, L~), a G (—4, 4) and L~ ^> L so that isometry along x~ direction 
holds approximately. Also we shall take a symmetric configuration around a = 0, hence ro = r(0) 
is an extremal point, r'(0) = 0. The induced metric reads as follows 

SVr = 2( _c t + c x) i SW = 4 + c 2 R r'(of . (2.3) 
From these expressions, the Nambu-Goto action takes the following form 

s = [ L/2 da (4 - 4) 1/2 (4 + 4 - V) 2 ) 1/2 • (2.4) 

y lixa Jo 

The energy is a first integral of motion, from which the following equation for the profile r(a) can 
be extracted 

f (A; 4(4-4) -1) , (2.5) 



r'(a) 2 



where k is an integration constant. Let us assume that the r.h.s. of ( |2.5| ) does not vanish at any 
location r G (r#,oo) with r# the location of the horizon, while, on the other hand, we assume 
cr(Th) — oo. Then, it necessarily holds that r = and the Wilson line extends symmetrically 
from r = oo down to r#. With these assumptions, which we must verify case by case, the profile 
can be obtained from 

Cr dr_ 

'r H c x (fc4(4-4)-l) 



= I ~ TT-T-^ TT7I • (2-6) 



In particular, the integration constant k is linked with L by the relation cr(oo) = -|. Going to a 
dimensionless radial coordinate p = r/rn, this is 



oc 



L = 2r H I ^ d ± -. (2.7) 

Ji cx (A:4(4-4) -1) 1/2 

The prescription in LRW for q calls for the leading behavior with L in the limit LT <C 1. This is 
clearly related to the limit k — > oo, i.e., 

2rH r CRdp + 0(k-^). (2.8) 



3 = ZeL f°° ^ ( Cx ~ c ^ Cx Cr dp (2 9) 

~ V2na> Ji (A;4(4-4)-l) 1/2 ' 

We must still subtract the contribution corresponding to the self-energy of the quarks. This is given 
by the Nambu-Goto action for a pair of Wilson lines that stretch straight from the boundary to 
the horizon, 

So = -^~, / ^(4-4)^ dp. (2.10) 



VkJi 4 (4 -4) 1/2 

Using (|2.5| ) and (p76|) , we can write the action as follows 



V2na' J i 
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To leading order in k 1 , taking into account that L is given by (|2.8|) , Sj = S — S reads ]T7| 



Si 



c R dp 



c x ( c x 



c 2 T yi 2 



(2.11) 



From here, let us extract an expression for the jet quenching parameter. We find it convenient to 
define 

<%{P) = 4-%{P) , c 2 x {p) = ^c 2 x {p) , c 2 R (p) = A R c 2 T (p) , (2.12) 

^r 

where the dimensionless quantity A R reads 



A; 



7-p 
r H 



1/2 



(2.13) 



A being the 't Hooft coupling in the p + 1 dimensional dual gauge theory. From these formulas, 
( |1.1| ) and (|1.2| ), we obtain 



1 



Q 



v 2rrA V •*') 



6— p 



cr dp 



c 2 T y/ 2 



(2.14) 



As it stands, this formula calls for a translation of r# in terms of the field theoretical quantities. 
In the case of non-rotating backgrounds we can provide a more explicit solution. For this class of 
metrics the Hawking temperature is given by the standard formula 



T 



4'(r) 



4tt v /4(^4( 



(2.15) 



r=r H 



Using this definition of the temperature, we can solve for ru/a' as follows 

4tt\/AT 



r JL 
a' 



c 2 T '(p) 



p=i. 



(2.16) 



and replace it in ( |2.14j ) to arrive at the formula 
1 



Q 



V2 



71 



167T 2 



9 -. 6 ~P 



c 2 T \l) 



T 2 (T 2 X)' 



cr dp 



i c 2 x {°\ 



c 2 T y/ 2 



(2.17) 



This expression is invariant under reparameterization of the radial coordinate (upon suitable change 
of the integration limits). Indeed, the dependence on T and A for generic p also coincides with the 
discussion in LRW. For example, as it stands, it can be directly used to extract the quenching 
parameter for the thermal deformation of Witten's D4-brane background. In this metric, the fifth 
dimension has been compactified to a circle of radius £. Hence, the four dimensional effective 
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coupling is A = X/£ = Ana>sYMN c . Therefore we may write for the effective quenching parameter 
the following expression 

g~ 14.26 cT 3 a SY MN c , (2.18) 

where c = IT is the ratio of radii of the thermal and Kaluza-Klein circles. The particular value c = 1 
signals the confinement/deconfinement transition temperature T — T c — l/£ ||25|| . Therefore, strictly 
speaking, ( |2.18| ) is valid for c > 1, or T > T c . The numerical factor in ( |2.18j ) is the result of ( |2.17| ) 



87rp 3 / 2 and 



3vr) 



-1/2 -3/2, 



1-p 



-3\ 



with p = 4 after inserting c 2 ^ = V S7rp 

]2"rj| . For standard values a s = ~ and iV c = 3, we get g = 0.87(2.93) GeV 2 /fm for c = 1 and 
T = 200 (300) MeV. These values are just slightly smaller than those in LRW. Yet, the 5d origin 
of ( |2.18|) is reflected in the linear dependence in the 't Hooft coupling. In the following sections we 
shall apply the expressions ( |2.14j) and ( |2.17| ) to another couple of relevant backgrounds. 

Another interesting use of the renormalized expression ( 2.17 ) is the possibility to extend the 
analysis to quark sources of finite mass. Indeed, from the point of view of the derivation, there 
is nothing peculiar about the integral upper limit being at the boundary, and it could equally 
well extend up to a finite value p m = r m jru < oo. The physics then is dual to a geometrical 
setting in which the ends of the fundamental string are attached to a probe brane that is placed 
in the above background at a fixed distance from the stack, r m = ma', set up by the mass of 
the quarks in the fundamental representation ]27 , 28[j . It is quite evident from the analytic form 



of the formula ( |2.17| ) that cutting the value of the integral will decrease the denominator, and 
hence enhance the value of q. Popular scenarios include cases where a probe D7 (alternatively a 
D6 or a D8) are placed inside D3 (respectively D4) backgrounds. Plotting q as a function of the 




Figure 1: Evolution of q with m for finite mass quarks. The solid curve represents the D3 background 
and the vertical axis is qj (\/asYM N c T 3 ) which asymptotes to 18.87... as in LRW. The dashed curve 
represents the D4 background. The vertical axis is now q/(casYM N C T S ) which asymptotes to 14.26. . . . 
The horizontal axis is p m = ma' /th- In each case, the quotient a' jru can be read off from ( 2.16j ). 



mass of the quarks, we observe a very weak dependence until the mass is rather low, where the 
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approximations are questionable. For example, for the case of the D3 background, the horizontal 
axis is p m = m/(2^l 2 ^a SY MN c T) ~ m/2.7 GeV for a SY M = 1/2, N c = 3 and T = 200 MeV. In 
order to increase q by a 10% we must lower p m until a value of ~ 10, hence m ~ 27 GeV. 



3. q at finite coupling 



The AdS/CFT correspondence is a statement that goes beyond the classical limit of string theory. 
In this limit, it maps classical solutions of supergravity to quantum field theory vacua in the strong 
coupling limit, A — > oo. Corrections in A -1 are in direct correspondence with those in powers of a' 
in the string theory side 4 . In this paper we shall use the solution given in |23|, 30 1 corresponding to 
the a' corrected near extremal D3-brane. The relevant pieces of information of this solution can be 
casted as follows 

C 2 T (P)=P 2 (1-P~ 4 )(l + 1T( P ) + ...), 
c 2 x (p)=P 2 (l+lX(p) + ...), 

clip) = p- 2 (l -P'THI + 7 R(p) + ...), 



to first order in 7 



CM f„J/r>2\3 



a'/R 2 ) 3 ~ 0.15 A 3 / 2 . Here, following fl2.13|) , we have already extracted 



Ar = R 2 jr 2 u factors, with R 2 = \/~\a' . The intervening functions read 

1225 o 695 



T(p) 
X(p) 
R{p) 



-75p- 4 
25 



16 



-P 



16 



■P 



-12 



16 

75p- 4 



p- 8 (i+p- 4 : 



(3.1) 



, 1175 o 4585 12 



16 ' 16 

Inserting this data in (|2.17 ) with p = 3 and expanding to first order in 7 gives 



q(\)=q(0) 1-7 



^-+45 
2a 



(3.2) 



with q(0) as given in [rB|, and 

dp 



r (5/4) 



^M^i r (3/4) ' 

R{p) - X{p) (p 4 + 2) + r(p) (p 4 -l) 

P 4 -1 



dp 



30725 7T 3 / 2 
924 v / 2~r(3/4) 2 



(3.3) 
(3.4) 



'Corrections in a' to timelike or spacelike Wilson lines used to compute the qq potential have appeared in 
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Evaluating yields 

q(X) =q(0) (1 - 1.7652 A~ 3/2 + ...) . (3.5) 

Therefore, we see that finite coupling corrections tend to diminish the value of the quenching 
parameter. For example, taking N c = 3 and olsym = 1/4, 1/2 and 1, we get a reduction factor of 
6, 2 and 0.8 percent respectively. Note that the decrease in the jet quenching parameter towards 
weak coupling is suggestive of a smooth interpolation between the strong coupling regime and the 
perturbative results. Obviously, the computation of higher order corrections would be necessary to 
put this conclusion on more solid grounds. 



4. q at finite chemical potential 

The near horizon metric of a rotating black D3-brane with maximal number of angular momenta 
reads as follows [[H], with the conventions of f3~3fl : 

/ 2 \ i 3 

ds 2 = VA y-H^fdt 2 + f~ l dr 2 + —dx-dxj+-j=Y^ R2R i [ dv 1 + "i(d<f>i + A'dtfl , (4.1) 



where V\ = cos#i, v 2 = sin^i cos #2? v z — sin^i sin 9 2 , and 7i = HiH 2 H 3 , where 



HI 



H 



1 + ^ 

j.2 ' 



f = — "H - E 



3 2 
Vt 



i=i 



(4.2) 



Upon Kaluza-Klein reduction, this becomes a charged AdS black hole solution of M = 2 U(1) 3 R 
supergravity, where A % plays the role of the gauge field. The holographic gauge theory corresponding 
to (3~T) is Af = 4 SU(N) supersymmetric Yang-Mills at finite temperature and with a chemical 
potential for the U(1) R symmetry. It will be convenient to trade the nonextremality parameter p 
for the horizon radius, r = r H given as the largest root of f{r H ) = 0, i.e., 



and define the adimensional quantities 



P = r f 2 n(r H ) 



4- , = 



2 ' 



As usual, we go to dimensionless variable p = r/r H and find 
Hi(p) = 1 + Kip~ 2 , 

Finally, the relevant functions entering the formula Q2.14 ) can be easily extracted from 



f(p) = ^ {p 2 n( P ) - P - 2 H(i)) = J-/( P ) . 

^r 



(4.3) 



(4.4) 



(4.5) 



c 2 T (p) 



A/ 



H 



1 3 



K; L Hi 



[H 



c 2 x (p) = Va~p 2 , c 2 R ( P ) 



f 



(4.6) 
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The factors in the metric depend on the internal angles. However, the terms above conspire to give 



c R dp 



1 



1 c xV c x 



H(l) k 



dp{p W) 



i 



-1/2 



(4.7) 



where all information about the internal angular coordinates has dissapeared. Now, given that the 
Hawking temperature of this solution is given by |31| 



(4.8) 



with R 2 = \/Aa', substituting in ( |2.14| ) we find the answer 



q(Ki 




2 + Ei=i - U i= i « 



-1/2" 



(4.9) 



In order to analyze this result, it must be recalled that the domain of thermodynamical stability 

is bounded by the inequality 5 n\ + k 2 + — Hi k 3 < 2. 



of this solution 



fsee also 



Plotting the right hand side of (|4.9|) numerically, we find the curves that are shown in figure ^. The 



V «SYM N c T 3 



60 
40 
20 



0.5 



1 



K\ =K2 

K\ =K2=K3 



1.5 



Figure 2: Plot of q as function of different combinations of charges Kj. The charges which are not varied 
are set equal to zero. The stability bound chops the lower curves, but leaves the upper one unbounded. 
At the origin the value is ~ 18.87 as found in LRW. 

jet quenching parameter raises its value for nonzero charges. The increase is not monotonous along 
the whole space of charges. In fact, though hardly noticeable, it changes the sign of the slope along 
the line Ki with ^2 = ^3 = 0. To see this better, we have zoomed the interval K\ e [1.4, 2] in figure 
|3], where we can see that the change in slope sign happens around Ki = 1.8. 

5 A wrong sign in this expression in a previous version of this paper has led to wrong plots that we have corrected 
in this version. 
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Figure 3: Zoom of the curve q{n\) showing the value of the turning and end points. 



A comparison with recent results in the literature is in order. We certainly agree and go beyond 
the perturbative analysis of Caceres and Guijosa Also, when restricted to the cases examined 
by Lin and Matsuo pOf (one charge) and Avramis and Sfetsos |2l] (one or two equal charges), we 
find qualitative agreement within the range of stability. We may also examine less symmetrical 
configurations on a 3 dimensional plot in figure 




Figure 4: Evolution of q in a sector of the plane (kj, K2 = K3). The curve is chopped by the thermody- 
namical stability bound. 



Besides analyzing the qualitative behaviour of q(ni) numerically, in order to compare with other 
approaches, it may be of interest to perform an expansion in powers of quantum field theoretical 
magnitudes. For this purpose, it is relevant to recall how the thermodynamical magnitudes are 
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related to geometrical quantities. In particular, the density of physical charge and chemical potential 
are given respectively by |35|, |37| 



i=l 

3 



^■Mu=^^n( i +*.) i/2 . <«■») 



where T = r^/(7ri? 2 ). From these expressions and ( |4.8| ), we should invert Kj in terms of pi and 
T for the canonical ensemble and in terms of //j and T for the grand canonical ensemble. This is 
dificult in the general case, so we may simplify for equal or vanishing values of For example, 
taking K\ = k and k 2 = K3 = the two (inverse) expansions yield 

K C = z-e+^e+- , (4.12) 



4 

^gc = C + C 2 + ^C 3 + -, (4.13) 



with 

2 



Expanding (|4.9| ) in powers of k and inserting these series we obtain for the canonical and grand- 
canonical ensemble respectively the following results 

Qc(p) =q(0) (1 + 0.63^ -1.08^ + 2.83^ + . ••) , (4.15) 
q GO {n) = q(0) (1 + 0.63 C + 0.18 ( 2 + 0.0Q( 3 + ...) . (4.16) 



This expansion fully agrees with the one in [21] upon rescaling £ — > 2£ and C — ► 2£ 



The computation of Wilson loops in thermal AdS backgrounds is a promising line of research, 
of which the q computation is a salient example. The identification of other observables within the 
AdS/CFT framework (see, for a recent example, ||38|| ) that can be confronted with experimental 
data is an urget challenge. It seems clear to us that there are several avenues for further exploration. 
Among these, the extension of our results to less supersymmetric backgrounds is of neat interest. 
We hope to report on these issues in the near future. 



6 The rescaling by 2 in £ is to be traced to the normalization of the gauge fields in |g5[, which is \/2 larger than 
usual. We thank Spyros D. Avramis for pointing this and some calculational errors in a previous version. 
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